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Abstract of Research Objectives 

The objective of this work is the comprehensive analysis of high- resolution 
atmospheric spectra recorded in the middle- infrared region to obtain simultaneous 
measurements of coupled parameters (gas concentrations of key trace constituents, 
total column amounts, pressure, and temperature) in the stratosphere and upper 
troposphere. Solar absorption spectra recorded at 0.002 and 0.02 cm 
resolutions with the University of Denver group's balloon-borne, aircraft-borne, 
and ground-based Interferometers and 0.005 to 0.01 cm resolution solar spectra 
from Kitt Peak are used in the analyses . 


Summary of Progress and Results 

The major recent focus of this work has been the analysis of the very high 
(0.002 to 0.003 cm -1 FWHM) resolution stratospheric solar spectra recorded during 
three recent balloon flights. The spectra reveal many new features of several 
important trace gases, such as HN0 3 and C0F 2 , as well as new details concerning 
the absorptions by Ci0N0 2 and H0 2 N0 2 . Comparison of the new spectra with line- 
by-line simulations shows that previous spectral line parameters are often 
inadequate, and that new analysis of high resolution laboratory and atmospheric 
spectra and improved theoretical calculations will be required for many bands. 
Preliminary versions of several sets of improved line parameters have been 
developed and compared with the balloon-borne spectra. A laboratory analysis 
of the v 6 band of C0F 2 has been completed. Isotopic ratios of 0 0 0 and 
16 0 is 0 16 0 relative to normal ozone have been measured based on recent improvements 
in ozone line parameters in the 10-^m region. Normalized to standard terrestrial 
isotopic ratios in ozone, both forms of heavy ozone are enriched a b°ve ^ ^ ^m 
altitude with about a factor of two larger enhancement in 16 0 16 0 18 0 than 0 0 0. 

In preparation for the Network for the Detection of Stratospheric Change, 
0.02 -cm -1 resolution solar spectra recorded at Mauna Loa have been analyzed. 
Simultaneous total columns have been determined for 13 atmospheric gases , and 
the daytime diurnal variations of the NO and N0 2 total columns have been 
measured. The results indicate that Mauna Loa is a favorable site for infrared 
monitoring of atmospheric gases. 

Analysis of the Kitt Peak ground-based solar spectra has continued. A series 
of 19 spectra recorded during the night using the Moon or during the day using 
the Sun as a source have been used to measure the daytime and nighttime variation 
of the N0 2 total column. A continuous increase is observed during the day 
(lifetime about 10.4 hours) followed by a rather rapid decrease during the night. 
At sunrise the N0 2 total column drops by about a factor of three. The 
measurements have been compared with values calculated with a 1-D diurnal 
photochemical model. Spectra recorded between December 1980 and May 1988 show 
that the CHCiF 2 (CFC-22) total column increased at an average annual exponential 
rate of 7.8% ± 1.0% (2cr). Compared with other atmospheric data, these 
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measurements indicate that CFC-22 is increasing at a more rapid rate than either 
CFC-11 or CFC-12, the two most abundant chlorof luorocarbons , but that the rate 
of CFC-22 increase is likely to have declined over the past few years. A 
systematic discrepancy between the absolute values of infrared and in situ CFC- 
22 measurements is noted. 
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A. Zonally-averaged Model of Dynamics, Chemistry and Radiation for the 

Atmosphere 

B. Principle Investigators : Nien-Dak Sze and Malcolm Ko 

Co- investigators : Arthur Y. Hou, Jose M. Rodriguez, Hans-R. Schneider 
and Debra K. Weisenstein 

Atmospheric and Environmental Research Inc., 840, Memorial Drive, 
Cambridge, MA 02139 

C. Abstract of Research Objective 

The overall goal of the project is to obtain an understanding of the 
processes that control the distribution and abundance of stratospheric ozone 
and its sensitivity to either natural or anthropogenic influences. The 
scientific objectives, which are designed around the two-dimensional zonal 
mean (2-D) modeling capability at AER, are : 

• To assess the limitation of 2-D models, to refine and enhance our current 
(2-D) modeling capabilities and understanding of the stratosphere through 
model simulations and validation of model results using observed data. 

• To develop an interactive 2-D model with interactive chemistry, dynamics 
and radiation. 

• To assess the susceptibility of stratospheric ozone to natural and 
anthropogenic perturbations. 

• To study the processes that give rise to the Antarctic Ozone Hole and 
assess their global impact. 


D. Summary of Progress and Results 
Ozone Budget 

We used result from the 2-D model to clarify the roles of chemistry and 
transport in determining the concentration of O 3 in the 2-D model (Ko et al, 
1989). We found that at the tropical lower stratosphere, the concentration of 
O3 is determined by a balance between chemical production and transport out of 
the region. At high latitudes, the balance is between local chemical removal 
and transport into the region. Our study help to explain the response of 
model ozone to chemical perturbation and clarify the limitation of 1-D model 
results . 

Assessment Calculations and Ozone Depletion Potentials 

The calculated O3 behavior in the future atmosphere from our assessment 
calculations were used in a study to examine the effect of redistribution of 
O3 on climate (Wang et al, 1989). We performed a number of assessment calcu- 
lations and provided inputs for the preparation of several reports including 
the NASA Reference Publication report 1208, "Present State of Knowledge of the 
Upper Atmosphere 1988 : An Assessment Report” and the UNEP/WMO report. We 
also used the 2-D model to calculate the Ozone Depletion Potential of a number 
of CFC substitutes and examined how the calculated values are affected by the 
transport parameters in the model. 

Interactive Model and Dynamics 

The interactive model is now operational. We have used the version with 
simplified chemistry to study the model sensitivity to Kyy (Schneider et al 
1989) and the feedback effects from change in O 3 . Our model simulations 
showed that the temperature and dynamics feedback from change in local heating 
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could compensate for about 1/6 of the calculated O 3 decrease. The version of 
the model with full chemistry has been used to simulate the present day atmos- 
phere and perform a number of steady state assessment calculations. The 
emphasis of the work is related to our theoretic study that clarify their role 
of Kyy in controlling the diabatic circulation (Hou and Ko, 1989). 

Antarctic Ozone 

Part of 1988 was spent in the analysis of the data from the Airborne 
Antarctic Ozone Experiment. Analysis of the data resulted in three publi- 
cations that discussed the effect on global 0 3 content from export of 0 3 -poor 
air from the vortex, the so-called dilution effect (Sze et al , 1989), the com- 
parison of simulated trace gas distributions in the vortex with the AAOE 
observations (Rodriguez et al, 1989) and the decadal trend of the ozone 
behavior in the polar vortex (Ko et al, 1989). 
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A. Acquisition and Utilization of Archived Satellite Data from the NASA 
Upper Atmospheric Research Program : Validation and Enhancement of a 
Two-dimensional Zonal-mean Model 


B. Principle Investigators : Malcolm K. Ko and Nien-Dak Sze 

Co- investigator : Arthur Y. Hou, Hans-R. Schneider and Debra K. 

Weisenstein 

Atmospheric and Environmental Research Inc., 840, Memorial Drive, 
Cambridge, MA 02139 


C. Abstract' of Research Objective 

The objective of the program is to acquire satellite data from the NASA 
Upper Atmosphere Research Program and to organise them in a form that is 
suitable for use in support of the two-dimensional modeling work at AER. 


D. Summary of Progress and Results 
Acquisition of Data 

We have acquired and organised five data sets. These are the BUV O3 
data; the SBUV 0 3 data; the TOMS 0 3 data; the LIMS data on temperature, 0 3 , 
NO 2 , H 2 O, HN0 3 and geopotential height; and the SAMS data on N 2 O and CH 4 . The 
data were interpolated into the AER model grid and processed that they are 
suitable for use as input for model studies and for comparison with model 
results . 

Calculation of Diabatic Circulation Using Observed Temperature and 0 3 

The zonal -mean 0 3 mixing ratio from SBUV was used together with the 
observed temperature from NMC to obtain the heating rate using the AER 
radiation code. A seasonal zonal-mean circulation was calculated from the 
diabatic heating rates. 

Interannual Variation of Ozone 

Using the 0 3 concentrations and observed monthly mean temperature from 
NMC for the four year period (1978-1982), radiative heating rates were com- 
puted for each month of the four years. The circulations diagnosed from the 
heating rates were used in the 2-D chemistry- transport model to calculate the 
0 3 behavior. The model calculated inter- annual fluctuations was compared with 
the observation. Our model calculation suggested that the observed inter- 
annual variability in 0 3 is a response to fluctuations in the transport 
circulation in the lower stratosphere. Furthermore, the year-to-year 
variations in the circulation are a results of the temperature fluctuations. 
Finally, the change in temperature cannot be explained by the fluctuations in 
0 3 and is more likely a results from the changes in eddy heat transport. 
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Utilization of N 2 O data from SAMS 

The atmospheric lifetime of N 2 0 is determined by the photochemical 
removal of the gas in the lower stratosphere. Previous calculations using 
balloon observations of N 2 0 as data and/or for validation obtained a lifetime 
for N 2 0 of about 150-175 years. Data for N 2 0 is available from the SAMS 
instrument on Nimbus -7 covering the region from 50° S to 70°N, 25 km to 50 km 
where 80% of the removal occurs. Using the observed N 2 0 data to check against 
the model calculated distribution, we concluded that calculated lifetimes as 
short as 114 years are consistent with the data. 

Photochemical Studies using LIMS data 

We incorporated the observed O 3 , H 2 0, and temperature from LIMS as input 
into our diurnal model. In this mode of operation, the values for local 
temperature, O 3 and H 2 0 are assigned the observed values and the 
concentrations for the other long-lived trace gases are taken from previous 
model simulations. The model is then used to calculate the concentrations of 
the short-lived species assuming local photochemical equilibrium. The model 
results are compared with the LIMS observed N0 2 and HNO 3 . 


Publication in Preparation 

Schneider, H-R, M.K.W. Ko, C.A. Peterson and E. Nash, Interannual Variation of 
Ozone : Interpretation of 4 years of Nimbus 7 SBUV Observations . 

Ko, Malcolm, N-D Sze and D. K. Weisenstein, Validation of Atmospheric 
Lifetimes for Source Gases Using N 2 0 Data from SAMS. 
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Studies of the antarctic ozone hole: 
Analysis of perturbations to global ozone trends 

Jose M. Rodriguez, Malcolm K. W. Ko, Nien Dak Sze 


Atmospheric and Environmental Research, Inc. 
840 Memorial Drive 
Cambridge, MA 02139 


Research Objectives 

The general objective of this program is to refine our understanding of 
the chemistry of the antarctic ozone hole, and to study the impact of polar 
phenomena (heterogeneous chemistry, dilution effect, denitrification) on the 
global ozone trends. The specific goals are: 

• Further studies of the dilution of the antarctic ozone loss using 
our two-dimensional chemistry transport model; comparison to other 
model results and observational data. 

• Incorporation of heterogeneous processes into our 2-D model, both in 
the polar regions and in the global atmosphere, in a manner 
consistent with available laboratory data and field measurements. 

• Assessment of the impact of heterogeneous reactions both for the 
present and future atmosphere; assessment of the future behavior of 
the antarctic ozone hole and its Impact on global ozone trends; 
assessment of the coupling among ozone and temperature 
perturbations, circulation changes, and dynamical redistribution. 


Progress and Results 

Our work during the first year has concentrated on three general areas: 
a) studies of the photochemistry of the Antarctic stratosphere; b) assessment 
of the Impact of high- latitude chemistry on the H0 X and O 3 budgets; c) impact 
of high- latitude and global heterogeneous chemistry on ozone trends. The main 
results of our research are as follows: 

• We have assessed the impact of new kinetic data for formation and 
photolysis of C ^202 on the partitioning of chlorine species and 
calculated ozone reductions within the Antarctic polar vortex. 
Calculations Indicate that the observed CiO and CiN 03 column 
densities near the edge of the vortex are consistent if we adopt the 
new kinetic data. Decreases in the observed column densities of HCi 
and C^N 03 towards the core of the vortex suggest an increase in C£0 
with latitude, with a correspondingly larger calculated reductions 
of ozone. The ozone reductions calculated with the new kinetic data 
are still consistent with observations, if we consider the spread in 
both the ozone and CiO measurements and mixing of air within the 
vortex . 

• We have incorporated heterogeneous chemistry in our 2-D model, both 
on Polar Stratospheric Clouds (PSCs) and on the global sulfate 
layer, 

• Model results indicate that the heterogeneous reaction of CiN 03 with 
H 2 O (ice, aerosol) could be the primary source of H 0 X at high 
latitudes winter (greater than 60°) if the atmosphere has been 
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depleted in N0 X but not in HN0 3 , and if the equivalent first-order 
rate of this reaction if about 10' 6 s . Inclusion of this reaction 
in our model Increased the calculated HO x abundances by almost a 
factor of 10, with corresponding changes in the partitioning between 
HCi and other chlorine species. Ozone loss below 22 km occurred in 
this case primarily through chlorine and HO x catalytic cycles. 

We have started the assessment of the impact of heterogeneous 
chemistry by considering two cases: 1) heterogeneous reactions 
occurring in a denitrified atmosphere on PSCs at high latitudes 
during winter; 2) heterogeneous reaction of HC,? with CiN0 3 occurring 
year-round in the global sulfate later, with a reaction efficiency 
of about 10 , Both cases yield calculated ozone trends for the 

past 20 years consistent with those derived by the NASA Ozone Trends 
Panel. However, the calculated future ozone trends are dramatically 
different. Possibility (1) yields ozone reductions in the year 2060 
of order 10% relative to 1985 at high latitudes, winter. The 
assumptions in (2) result in calculated ozone reductions of order 20 
- 30% for the same year, season, and latitude. Uncertainties in the 
kinetic data for heterogeneous reactions indicates the need for 
further sensitivity studies, and for validation of model results by 
comparing to data on other trace species (ie., CJIO) . 


Conference Presentations 

Rodriguez, J. M. , M. K. W. Ko, and N. D. Sze (1989) "Impact of new kinetic 

data for the chemistry of the Antarctic stratosphere," presented at the 
Spring Meeting of the American Geophysical Union, Baltimore, MD, May 7- 
12, 1989. EOS, 70, 299. 

Rodriguez, J. M. , M. K. W. Ko, N. D. Sze, and D. Weisenstein (1989): "Global 
implications of heterogeneous chemistry: Assessment utilizing a two- 
dimensional model," presented at the Fifth Scientific Assembly of the 
International Association of Meteorology and Atmospheric Physics (IAMAP), 
Reading, England, July 31 -August 11, 1989. 


Publications 

D. W. Fahey, D. M. Murphy, K. K. Kelly, M. K. W. Ko, M. H. Proffitt, C. S. 

Eubank, G. V. Ferry, M. Loewenstein and K. R. Chan (1989) "Measurements 
of nitric oxide and total reactive nitrogen in the Antarctic 
stratosphere: Observations and chemical implications," J. Geophvs . Res., 
in press. 


Publications in preparation 

Rodriguez, J. M. , M. K. W. Ko and N. D. Sze (1989): "Partitioning of chlorine 
species inside the Antarctic polar vortex." 

Rodriguez, J. M. , M. K. W. Ko, N. D. Sze and D. Weisenstein (1989): "Global 
implications of heterogeneous chemistry." 
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Research Summary 1988-1989 


A. Title of Research Task 

Interactions of Atmospheric Chemical, Radiative, and Dynamical Processes: Research 
Studies with and Further Development of the LLNL Two-Dimensional Global Atmospheric 
Model 


B. Investigators and Institutions 

Principal Investigator: Dr. Donald J. Wuebbles 

Co-Investigators: Dr. Peter S. Connell 

Dr. Keith E. Grant 
Dr. Douglas E. Kinnison 


C. Research Objectives 

The primary goals of this project are the further development of the LLNL two- 
dimensional chemical-radiative-transport model of the troposphere and stratosphere, and 
the application of this model to research studies on the middle atmosphere. Model 
capabilities are being extended towards improving understanding of the key processes that 
control the distribution of ozone and other important species in the global atmosphere. 
Research studies are primarily oriented at assessing the past and possible future changes in 
ozone from natural and human-induced causes, and to improve understanding of 
observations of global trace constituent distributions. 


D. Progress and Results 

The two-dimensional model has been applied to a number of research studies during 
1988 and 1989. LLNL scientists have also participated in several international assessment 
efforts in support of NASA during this time period. Past trends in ozone and temperature 
were analyzed with the model for the Trends Panel Report (we also contributed to the 
writing of this report being published through NASA and other organizations). For the 
recent Scientific Assessment of Stratospheric Ozone: 1989 (being published by WMO), we 
not only provided model calculations for trace gas scenario effects on ozone and an analysis 
of ozone depletion potentials for CFCs and other compounds, but also coordinated the 
preparation and writing of one of the four chapters. In 1988, we helped coordinate and 
participated in an international two-dimensional intercomparison of atmospheric models 
sponsored by NASA. 

Research studies with the two-dimensional model include: an analysis of trends in 
ozone and temperature over the past thirty years, including a careful evaluation of the 
possible effects on ozone resulting from atmospheric nuclear testing in the early 1960s; an 
analysis and sensitivity study of potential effects of high speed commercial aircraft on 
stratospheric ozone; an analysis of carbon- 14 and other tracers produced from nuclear tests 
to determine their usefulness for validating the treatment of troposphere-stratosphere 
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exchange processes in the 2-d model; and an analysis of the 1 1-year solar cycle for its 
effects on upper stratospheric ozone and temperature. 

The two-dimensional model has been continued to be improved and refined. 
Additional layers have been added to the model to extend it into the lower mesosphere and 
to improve the resolution in the troposphere. The solar and infrared radiative submodels 
have extensively revised to provide improved treatment of scattering and cloud-related 
processes. Determination of upper stratospheric cooling rates has also been improved. At 
the same time, computational speed has been improved, such that the current model runs as 
fast as earlier versions. 
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A. Two-Dimensional Model Studies of the Middle Atmosphere 

B. R. R. Garcia, National Center for Atmospheric Research, Boulder, CO 
S. Solomon, NOAA / Aeronomy Laboratory, Boulder, CO 

J. T. Kiehl, National Center for Atmospheric Research, Boulder, CO 

C. The two-dimensional model of Garcia and Solomon will be used to study several 
aspects of radiative-chemical-dynamical coupling in the middle atmosphere. The 
proposal consists of three components: Improvements and extensions of our present 
2D model; studies of radiative/photochemical coupling in the contemporary middle 
atmosphere; and studies of possible interactions among radiation, dynamics and 
photochemistry that may be important in the future atmosphere due to changes in 
CFCs, carbon dioxide, methane, and nitrous oxide. 

D. The 2D model has been improved by the addition of an IR transfer model developed by 
J. Kiehl. This model uses a band formulation to explicitly compute IR heating/cooling 
rates. The model domain had been extended to include die troposphere. Model 
dynamics have been reformulated to include explicit solution of the zonal mean 
momentum equation. The new dynamical formulation should allow studies of tropical 
dynamics not possible with the previous version wherein the zonal wind was assumed 
to be in thermal wind balance with the temperature field. A single-wavenumber 
planetary wave model has been coupled to the 2D model to study seasonal and 
interhemispheric differences in transport due to planetary wave activity. 

E. Papers on the response of the stratosphere to multiple anthropogenic perturbations 

(CFCs + C02 + CH 4 ) and on the new model dynamics are in preparation. 
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NASA RESEARCH SUMMARY 
Grant No. NAGW-1605 


A. Title of Research Task: 

Zonally- Averaged Model of Dynamics, Chemistry and Radiation for the 
Atmosphere 

B . Principle Investigator: 

Ka-Kit Tung, Professor of Applied Mathematics, University of Washington 

C . Abstract of Research Objectives: 

To construct a physically based 2-D model with coupled dynamics, chemistry and 
radiation; to use such a model to simulate the present date distribution of trace gases in the 
stratosphere and to assess future change. 

D. Summary of Progress and Results: 

We have completed the implementation of the diagnostic model using NMC 
temperature as input from which both the advective and diffusive transports are deduced. 
We have recently added a narrow-band radiative transfer model, which has an accuracy of 
better than 10% when compared to the line-by-line calculations. We have also completed 
the implementation of the chemistry module, with O x , HO x , NO x and CIO x cycles. 
The simulation of present date distribution of many species, including ozone, using the 
coupled model showed great improvements over previous generations of models. 

E. Publications: 

K.K. Tung and H. Yang: Dynamical component of seasonal and year-to-year 
changes in Antarctic and global ozone. J. of Geophys. Research, 22, 
12537-12559, 1988. 

K.K. Tung and H. Yang: Dvnamic variability of column ozone, J. of Geophys. 
Research, £2, 11123-11128, 1988. 

H. Yang, K.K. Tung and E.P. Olaguar: Two-D model simulation of ozone 
climatology and year-to-year variations, in Proceedings of Quadrennial 
Ozone Symposium, 1988. 

H. Yang, K.K. Tung and E.P. Olaguar: Nongeostrophic theory of zonally 
averaged circulation. Part II: E-P fluxes and isentropic mixing coefficients, 
J. Atmos. Sci., accepted 1989. 
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